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The dichotomy of the pion as QCD's Goldstone mode and a bound state of massive constituents 
is easily understood using the Dyson-Schwinger equations. That provides the foundation for an 
efficacious phenomenology, which correlates the pion's charge radius and electromagnetic form factor 
with its valence quark distribution function; and simultaneously provides a Poincare covariant 
description of the nucleon, its form factors and, more recently, meson photoproduction processes. 
This well-constrained framework can also be used to eliminate candidates for an extension of the 
Standard Model by providing the relation between current-quark electric dipole moments and that of 
' the neutron. 

^C) ' — A summary of two presentations, one by each author. 
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1 Introduction 

> 

A focus of contemporary studies in QCD is the development of an intuitive understanding of the 
spectrum and interactions of hadrons in terms of QCD's elementary excitations; i.e., quarks and 
gluons. Progress can be made by applying a single framework to the calculation of many observ- 
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ables. This facilitates a verification of necessary model assumptions, and the identification of robust 
correlations between a theory's keystones and hadron properties. Non-hadronic electroweak inter- 
actions provide the obvious test bed for any such approach because the electroweak probes are well 
understood and hence a given experiment yields immediate access to properties of the hadron tar- 
get. Thus constrained the framework can be used reliably to make predictions for other phenomena, 
even those far removed from the domain on which it was constrained; e.g., the properties of QCD 
at nonzero temperature and baryon density. 

Herein we supply a brief description of recent progress with the Dyson-Schwinger equations 
(DSEs) JH in these applications. A familiar DSE is the gap equation that describes pairing and 
condensation in low temperature superconductors; another is the Bethe-Salpeter equation, whose 
solution provides the mass and "wave function" of a bound state in quantum field theory. These 
are subjects in which the nonperturbative character of DSEs is paramount, to which we will return. 
However, at their simplest, the DSEs are a generating tool for perturbation theory: the weak 
coupling expansion of any particular DSE yields all the well-known Feynman diagrams. This is 
of immense help in studying QCD because it means there is little or no model dependence in 
the ultraviolet behaviour of calculated quantities. The key model-dependence is limited to the 
infrared domain, a property that can be exploited to probe the dynamics underlying confinement 
and dynamical chiral symmetry breaking (DCSB), which are QCD's signature nonperturbative 
phenomena. The recent successes and current challenges in this application are documented in 
Refs. gg. 

2 Dynamical Chiral Symmetry Breaking 

While the dynamical breaking of chiral symmetry in QCD is fundamental to the success of chiral 
perturbation theory: it is only owing to DCSB that m n = in the chiral limit and the scale of 
f n is set by the constituent-quark mass, understanding the origin of this phenomenon is outside 
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Figure 1: QCD Gap Equation or Dyson- 
Schwinger equation for the quark self energy, 
S. Here: D is the dressed-gluon propaga- 
tor; r, the dressed-quark-gluon vertex; and 
S = 1 / [ij • p + m + the dressed-quark 

propagator. 



the scope of the theory. This is where QCD's gap equation, depicted in Fig. |], finds immediate 
application 0. The solution provides the dressed-quark propagator in terms of the self-energy, X: 



S(p) 



l/[*7 • p + m + E(p)] , E(p) = 17 ■ p [^(p 2 ) - 1] + B(p 2 



(1) 



where m is the current-quark mass, and the functions A(p 2 ), B(p 2 ) are completely determined by 
the nature of the force between quarks. 

Employing a weak coupling expansion yields the perturbative results: A(p 2 ) ~ 1 and 
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where the ellipsis indicates that, like the zeroth and first order contributions, every higher-order 
term in the perturbative evaluation of B(p 2 ) is proportional to the current-quark mass. Hence at 
every finite order in perturbation theory the scalar piece of the quark's self-energy, B(p 2 ), vanishes 
in the chiral limit, m = 0. Therefore dynamical mass generation, and hence DCSB, is impossible 
in perturbation theory. 

Now the essentially nonperturbative character of DSEs becomes important. The self-energy 
appears in the denominator of the r.h.s. in Fig. p] but in the numerator on the l.h.s. This makes 
the equation nonlinear and hence its self-consistent solution can exhibit properties inaccessible in 
perturbation theory. Using simple models of the equation's kernel one can unambiguously establish 
that B{p 2 ) ^ is the favoured solution in the chiral limit if, and only if, the integrand provides 
sufficient support on the domain k 2 E [0,2] GeV 2 ; i.e., if the effective coupling in the infrared is 
strong enough A realistic, one-parameter model of the effective interaction || yields the mass 

function depicted in Fig. |2| It is important to remember that the existence of a nonzero solution 
in the chiral limit is a purely nonperturbative effect. Hence the domain on which the chiral limit 
solution and the u-quark solution are nearly indistinguishable is that on which nonperturbative 
dynamics is dominant in QCD: where they separate marks the beginning of the perturbative domain 
and this point is characterised by a length-scale of ~ 0.15 fm. 

The quark condensate is a fundamental fitting parameter in chiral perturbation theory. In DSE 
studies it is a calculated quantity that can simply be read-off from the large-p 2 behaviour of the 
dressed-quark mass function. In Landau gauge 
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where 7 m = 12/(33 — 2Nf), Nf = 4, is the anomalous mass dimension, and (qq)° is the gauge- 
invariant and renormalisation-point-independent vacuum quark condensate, which is easily evolved 
to the mass-scale ~ 1 GeV relevant to chiral perturbation theory. The one-parameter model of 
Ref. § yields 
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Figure 2: Quark mass function: AI(p 2 ) = 
B (p 2 ) I A(p 2 ) , in the chiral limit, and for 
the it-quark (isospin symmetry is assumed): 
m iGeV _ 5.5 MeV, and the s-quark: 
m s GeV = 130 MeV. The intersection of the 
line p 2 = M 2 (p 2 ) with a given curve gives 
the Euclidean constituent-quark mass ||. 



a value consistent with recent lattice simulations m. The condensate's value is a QCD analogue 
of the Cooper pair density in a BCS superconductor, and the value in Eq. (||) corresponds to 
a density of 1.8 fm~ 3 . Each sphere in a close-packed sea with this density would have a radius 
r/qg\ = 0.76 fm, which is just 15% larger than the pion's charge radius and 13% less that the 
proton's. This comparison underscores the importance of a proper description of DCSB in any 
attempt to explain low-energy phenomena. 

3 Mesons 

An aim of contemporary experiments is to explore the transition between the nonperturbative and 
perturbative domains in QCD. That requires momentum transfers at which a theoretical description 
of the reactions must be Poincare covariant. For mesons this means employing an homogeneous 
Bethe-Salpeter equation (BSE) to calculate the mass, and the amplitude that will be used in 
evaluating necessary matrix elements. (While QCD's gap equation is a DSE for the dressed-quark 
2-point function, a propagator, the BSE is a DSE for a 3-point function; i.e., a vertex.) There 
is a direct connection between the Bethe-Salpeter bound state equation and the vertices that 
appear in the Ward-Takahashi identities that are a true field theoretical representation of current 
conservation. The simplest of these identities relate the vertices to the dressed-quark propagator 
and therefore their fulfillment is only possible if there is a tight connection between the kernel in 
QCD's gap equation and that in the Bethe-Salpeter equations. 

It is a feature of quantum field theory that the Bethe-Salpeter kernel cannot be written in a 
closed form and hence all concrete calculations must employ a truncation. This is also true of the 
gap equation's kernel. It is commonplace to find calculations that ignore the constraints applied by 
the Ward-Takahashi identities, and employ kernels in the the gap equation and the BSE that are 
incompatible. Such studies necessarily violate the chiral symmetry constraints which the success of 
chiral perturbation theory has shown to be so important in low-energy QCD. It need not be thus. 

3.1 A Mass Formula 

There is at least one practical, systematic, symmetry preserving truncation of the DSEs II. It has 
been used |]] to prove Goldstone's theorem in QCD, to obtain quark-level Goldberger-Treiman 
relations that relate the scalar functions in the pion's Bethe-Salpeter amplitude to those in the 
dressed-quark propagator, and to derive a mass-formula for flavour nonsinglet pseudoscalar mesons: 

f 2 H m 2 H = -(q q )?M?. (5) 
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In this equation: M.^ = rnj} + m^ 2 is the sum of the current-quark masses of the meson's con- 
stituents, with C the renormalisation point; the electroweak decay constant is obtained from 

f H P, = Z 2 (C,A) / A (0 tr [[lT H ) t % ^S(q + )T H (q ] P)S(q^] , (6) 

with Z2(C,A) the quark wavefunction renormalisation, A the mass-scale in a translationally in- 
variant regularisation of the integral, q± = q ± P/2 (P^ is the meson's momentum and F H is its 
bound state amplitude), and, e.g., T 7r+ = (A 1 + iA 2 )/2, where {A J , j = 1, . . . ,8} are the Gell-Mann 
matrices, and (•)* denotes matrix transpose; and the in-hadron condensate is 



i<W>f = fHZ^K) / A (0I tr [(lT H ) t l5 S(q + )r H (q ] P)S(q^ . (7) 



The Z2 on the r.h.s. of Eq. (|6|) ensures that fjj is gauge invariant, and cutoff and renormalisation- 
point independent; i.e., that it is observable. (Equation. (||) is the field theoretical expression for the 
pseudovector projection of the pion's wavefunction at the origin in configuration space.) Further- 
more, the quark-level Goldberger-Treiman relations and Eq. (@) make plain that in the presence of 
DCSB the magnitude of the pion's leptonic decay constant is set by the constituent-quark mass. 
The same is true for (qq)^ , which, in the chiral limit, is identical to the vacuum quark conden- 
sate ||. The factor Z4 ensures that the in-hadron condensate is gauge and cutoff independent, 
and that its renormalisation-point dependence is precisely that required to ensure (qq)f -M( is 
renormalisation-point in-dependent. 

Equation (|B|) is written suggestively: it has the appearance of the Gell-Mann-Oakes-Renner 
relation, and it can be shown Q that for small current-quark masses it does indeed coincide with 
that formula. The new aspect of the equation is that it is valid independent of the current-quark 
mass of the constituents: the DSE derivation assumes nothing about the size of m^ 1 ' 92 . It has 
consequently been used to prove that the mass of a heavy pseudoscalar meson rises linearly with 
the mass of its heaviest constituent [10]. Equation (||) is a single mass formula that unifies the light- 



and heavy-quark domains. It also provides [11] an understanding of recent lattice-QCD data on the 
current-quark mass-dependence of pseudoscalar meson masses. 

The truncation scheme of Ref. M is the foundation for a phenomenological model that has 



been used to very good effect in describing light-quark mesons |g, 12]. That model is the only one 



to predict [13] a behaviour for the pion's electromagnetic form factor that agrees with the results 
of a recent Hall C experiment [14]. The large-Q 2 behaviour of the form factor can be obtained 



algebraically and one finds [|T5| Q F n (Q ) = const., up to logarithmic corrections, in agreement with 
the perturbative-QCD expectation. This result relies on the presence of pseudovector components 
in the pion's Bethe-Salpeter amplitude, which is guaranteed by the quark-level Goldberger-Treiman 
relations proved in Ref. Q. 

3.2 Pion's Valence Quark Distribution Function 

The DSEs provide a chiral-symmetry preserving, dynamical approach to QCD, which easily cap- 
tures the dichotomous nature of the pion as: 1) QCD's Goldstone mode; and 2) a bound state of 
quarks with large constituent-masses, and unifies the low- and high-Q 2 domains. They are therefore 
well-suited to a study of the pion's valence-quark distribution function, Uy(x). (This is a measur- 
able expression of the pion's quark-gluon substructure but it cannot be calculated in perturbation 
theory.) Since it targets are scarce other means must be employed to measure u v (x) and one ap- 
proach is to infer it from ttN Drell-Yan [p|, which yields Uy{x) oc (1 — x) for x ~ 1. However, a 
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DSE calculation predicts [17|: 



u v {x) ^(l-x) 2 , (8) 
at a resolving scale of /x 2 ~ 1 GeV 2 , and while this conflicts with the extant experimental result, as 



apparent in Fig. ||, it is consistent with the perturbative-QCD expectation [18]. The discrepancy is 
very disturbing because a verification of the experimental result would present a profound threat to 
QCD, even challenging the assumed vector-exchange nature of the interaction. The DSE study p7| 
has refocused attention on this disagreement, and is the catalyst for a resurgence of interest in 
Uy(x) and proposals for its remeasurement [19]. 



4 Baryons 

A direct analogy to treating mesons via the Bethe-Salpeter equation is to describe baryons using a 
Poincare covariant Fadde'ev equation. Of course, that equation also involves a kernel about which 
assumptions must be made in order to arrive at a tractable problem. Here the truncation scheme 
of Ref. H provides, a posteriori, a basis for treating baryons as quark-diquark composites using a 
Fadde'ev equation of the type proposed in Ref. pC| ]. Therein two quarks are always correlated as a 
colour-antitriplet diquark quasiparticle (because ladder-like gluon exchange is attractive in the 3 C 
quark-quark scattering channel) and binding in the nucleoli is effected by the iterated exchange of 
roles between the dormant and diquark-participant quarks. A first numerical study of this Fadde'ev 
equation was reported in Ref. |21], and following that there have been numerous more extensive 



analyses and applications; e.g., Refs. [22, which are reviewed in Ref. [24] 
4.1 Fadde'ev Amplitude Ansatz 

In developing an efficacious phenomenology it is possible to bypass solving the Fadde'ev equation, 
which can be a numerically intensive process, and employ a product Ansatz for the nucleon's bound 
state amplitude; an approach kindred to that which is still employed fruitfully in the study of meson 
properties, e.g., Ref. |plf |. The simplest Ansatz retains only a scalar (0 + ) diquark correlation and 
models the nucleon's amplitude as |2^j : 

V>(Pi,ai,7i) oc e clC2C3 [r 0+ (I(pi -P2);K)]ZZ A 0+ (K) [S(l;P)u{P)]* , (9) 
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Obs. 


Calc. 


{r v Y (fm 2 ) 


(0.87) 2 


(0.78) 2 


(r n ) 2 (fm 2 ) 


-(0.34) 2 


-(0.40) 2 


t^p Vr^JV / 


2.79 


2.85 


Ur> (Uiv) 


-1.91 


-1.61 


f^TL / f^p 


-0.68 


-0.57 


9wNN 


13.4 


13.9 




(0.93-1.06) 2 


(0.63) 2 


9A 


1.26 


0.98 


<r 2 ) (fm 2 ) 


(0.68 ±0.12) 2 


(0.83) 2 


9pNN 


6.4 


5.61 


9ujNN 


7-10.5 


10.0 



Table 1: Calculated values of a range of 
physical observables obtained [26] using the 
Fadde'ev Ansatz parameters in Eq. (12). 
The "Obs." column reports experimental 
values 0] or values employed in a typical 
meson exchange model |2£ 



where: (pi,cti,Ti) are the momenta, spin and isospin labels of the quarks comprising the nucleon; 
e ClC2 c 3 is the Levi-Civita symbol that gives the colour singlet factor; P = p\ + P2 + P3, K = p\ +P2 
and t = (—pi — P2 + 2p3)/3; A (K) is a pseudoparticle propagator for the scalar diquark formed 
from quarks 1 and 2, and r 0+ is a Bethe-Salpeter-like amplitude describing their relative momentum 
correlation; S, a 4 x 4 Dirac matrix, describes the relative quark-diquark momentum correlation; 
and u(P) is a free-nucleon spinor. The unknown functions are parametrised as (3~{x) = (1— e~ x )/x): 

A 0+ (K) = -Lf(K 2 /u; 2 0+ ), T 0+ (k;K) = J ^ T C ll5l r 2 ^(k 2 /u; 2 0+ ), (10) 
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s(£;P) = —r^/uA) i D - — ii 7 .e-e-pi D ) , (ii) 



where: C = 7274 is the charge conjugation matrix; is a calculated normalisation factor that 

guarantees an electric charge of 1/3 for the scalar diquark; is a similar factor that ensures 
the proton has unit charge; and R = 0.5 measures the ratio of lower to upper component in the 



nucleon's spinor in the rest frame, which is taken from Fadde'ev equation studies [11|. The model 
has three parameters: mj, the diquark's mass; and loq+ and u^p. They have physical intepretations: 
d + = l/m + is the distance over which a scalar diquark correlation can propagate inside the 
nucleon; l + = 1/uJq+ is a measure of the mean separation between the quarks in the scalar diquark; 
and = measures the mean separation between the dormant quark and the diquark. 

In the many applications of this Ansatz the parameters are typically determined by requiring a 



m + 


u + 




l/m + 


iM)+ 


1/Uty 


0.62 


0.79 


0.23 


0.32 


0.25 


0.86 



least-squares fit to the proton's electromagnetic form factor, as described in Ref. [26|. The procedure 
yields (in GeV or fm, as appropriate) 

(12) 

and it is plain that these values provide an internally consistent picture: l + < L,, which means 
that the dormant quark is spatially separated from the diquark; and d + < r p , the proton's charge 
radius, so that the diquark doesn't propagate outside the nucleon. The subsequent calculated results 
are then predictions and tests of the model's fidelity, which can be gauged from Table [|. 

4.2 Meson Photoproduction 

The same Ansatz is being used to study photoproduction of mesons from the nucleon. These pro- 
cesses are important for developing an understanding of the structure of nucleon resonances and in 
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Figure 4: t-channel 7r-exchange contribution 
to the 7r-photoproduction amplitude. In me- 
son exchange models the 77t*7t vertex is usu- 
ally considered momentum-independent. 



searching for "missing" resonances; i.e., those states predicted by constituent quark models that are 
hitherto unobserved. The aim of these calculations is to provide and constrain the input to meson 
exchange models, which are the tool that makes possible a comparison with data by subsequently 
incorporating details of the reaction mechanism. First results for to photoproduction were presented 
in Ref. and, while only the i-channel 7r-exchange mechanism was described, the results at for- 
ward angles, where this process must dominate, demonstrate the approach's promise. A calculation 
of the contributions to the cross-section from the s- and it-channel processes is almost complete. 

Meanwhile the analogous contribution to pion photoproduction, illustrated in Fig. |], has been 
calculated and the parameter-free results are depicted in Fig. |H[ The solid curve is the Born ap- 
proximation calculated using an efficacious meson exchange model |H|]. That model neglects the 
necessary momentum dependence of the 77r*7r vertex, for which the DSEs make a prediction. Includ- 
ing the calculated vertex yields the short-dashed curve in Fig. ||, which is not materially different. 
A posteriori, this justifies the meson exchange model expedient of neglecting quark-gluon sub- 
structure for these vertices. Of course, with the new DSE predictions, quantitative improvements 
are now possible. The dashed and dot-dashed curves in Fig. |5| were obtained using the calculated 
dipole width for the irNN vertex: F^nn^-^), which is larger (< 2) than that fitted to data in 



Ref. [31]. The effect is significant in this case. However, there are preliminary indications that the 
particularly soft form factor used in the meson exchange model is implicitly also accounting for off- 
nucleon-mass-shell suppression in the form factor. It is a known property (see, e.g., Ref. |25|) that 
vertices describing the interaction of three composite objects provide suppression when any one of 
the attached legs is off-shell. The DSE prediction for the strength of this effect will facilitate the 
qualitative improvement of meson exchange models by making possible the explicit representation 
of this phenomenon. 

4.3 Neutron's Electric Dipole Moment 

A well-founded description of the nucleon makes many applications possible, even constraining 
extensions of the Standard Model, which is an important focus of contemporary nuclear and particle 
physics. It has long been known that the possession of an electric dipole moment (EDM) by a spin- 
1/2 particle would signal the violation of time-reversal invariance. Any such effect is likely small, 
given the observed magnitude of CP and T violation in the neutral kaon system, and this makes 
neutral particles the obvious subject for experiments: the existence of an electric monopole charge 
would overwhelm most signals of the dipole strength. It is therefore natural to focus on the neutron, 
which is the simplest spin-i neutral system in nature. Attempts to measure the neutron's EDM, 
d n , have a long history and currently |$^] 

\d n \ < 6.3 x lCT 26 ecm (90%C.L.). (13) 

(NB. e/(2M n ) = 1.0 x 10~ 14 ecm. Therefore, writing d n = eh n /(2M), where M is the neutron's 
mass and h n is its "gyroelectric ratio," Eq. ([[^) corresponds to \h n \ < 6.0 x 10~ 12 .) 
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Figure 5: ^yN 



irN differential cross- 



section at Ej = 0.34 GeV as obtained solely 
from the Born diagram in Fig. |j. The solid 



line is the Born cross-section obtained 3C] 
using the meson exchange model of Ref. 1 31 ] 
whose parameters were fitted to data in a full 
T-matrix calculation. The remaining curves, 
discussed in the text, are calculated using 
DSE input. 



This experimental constraint on d n has been very effective in ruling out candidates for theories 
that enlarge the Standard Model. That is true because in the Standard Model the first nonzero 
contribution to a free quark's EDM appears at third order and involves a gluon radiative correction 
(i.e., 0(a s G%), for the same reason that flavour-changing neutral currents are suppressed: the GIM 
mechanism) so that |26|| 

d SM < 10 "34 ecm _ (M) 

This is seven orders-of-magnitude less than the experimental upper bound. However, the Standard 
Model is peculiar in this regard and candidates for its extension typically contain many more 
possibilities for CP and T violation, which are not a priori constrained to be small. Hence Eq. ( |i~3|) 
is an important and direct constraint on these extensions because Eq. ( |i"4| ) indicates that the 
Standard Model contribution to d n cannot interfere at a level that could currently cause confusion. 

Extensions of the Standard Model are typically used to predict current-quark EDMs, and to 
proceed from these to a result for d n one must have an understanding of the relation between 
current-quarks and constituent-quarks, and a reliable model of the neutron. The DSEs give both: 
the first from the gap equation, and the second from the Fadde'ev equation studies and their 



phenomenological application. Putting these elements together one arrives [26] at a quantitative 
relation between the current-quark's gyroelectric ratio and that of the neutron: h n ~ —91 hd- Here 
the very large magnifying factor: 91, owes its appearance to DCSB, which turns the current-quark 



into the constituent-quark, and it means that Eq. (13) applies the following bound: 

\h d \ < 7.4 x 1(T 14 . (15) 
(NB. Most extensions of the Standard Model predict h u <C hd-) This DCSB-tightened bound very 



much threatens the viability of a popular three-Higgs-boson model of spontaneous CP violation [33]. 
Furthermore, and importantly, the result in Eq. (|i~5| ) is independent of the model used to calculate 
hd and hence can be applied directly to constrain any extension of the Standard Model. 

5 Epilogue 

Owing to dynamical chiral symmetry breaking in QCD the pion appears as both a Goldstone mode 
and a bound state of massive constituents. The existence of at least one systematic, symmetry 
preserving truncation scheme makes a detailed understanding and explanation of this phenomenon 
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possible using the Dyson-Schwinger equations (DSEs). This is the starting point for a successful 
phenomenology of strong interaction phenomena, which quantitatively and qualitatively unifies 
the low- and high-energy regimes. Since the DSEs maintain contact with perturbation theory, 
the model-dependence is restricted to a statement about the infrared behaviour of the quark- 
quark interaction; i.e., the unknown nature of the long-range force in QCD. This remaining model- 
dependence is a virtue because it makes possible the correlation of observables via a parametrisation 
of this infrared behaviour and hence the use of experiments as a probe of the confining force. 

Acknowledgments: This work was supported by the US Department of Energy, Nuclear Physics Division, 
under contract no. W-31-109-ENG-38. 
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